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Abstract. In the beginning of this note we review the fundamentals of Grass- 
mann and flag bundles, stressing their key role in the geometrical formalism 
of nonlinear PDEs. Then we observe that the equation of involutive planes 
of the Cartan distribution — a nonlinear PDE naturally associated with any 
jet space — is equipped with a polar distribution. The equation plays the role 
of "empty equation" in the context of singularity equations, while the polar 
distribution is a source of contact invariants. In the last part, some examples 
of polar distributions are described explicitly. 
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Introduction 

The equation of involutive fc-planes of the Cartan distribution (Xfc(C), in this 
note) can be seen as a coordinate-free way of defining higher-order jet spaces out 
of the lower-order ones; as such, the idea behind it is very old, and it is impossible 
to point an adequate reference. Most of the textbooks [5J [TJ [2DJ [THl [B5] on the 
geometry of jet spaces and nonlinear PDEs, in one form or another, exploit it; it 
can be also found in the neighboring area of geometric mechanics disguised as a 
"constraint" in the sentence higher-order mechanics can be seen as a constrained 
first-order ones. However, it has never given a special attention as an equation 
per se, since its role is not fundamental for the study of the geometry of nonlinear 
PDEs and their (regular) solutions. 

Nonetheless, if singular solutions are brought into play, then the role of Ifc(C), 
where k is less than the number of independent variables, becomes not so choreo- 
graphic — a good reason, in the authors' opinion, for studying it in detail (Section 
[3u . Besides the context of singular solutions, whose main contributions come from 
Vinogradov, Krishchenko and Lychagin |BJ 00[IDl[nj [E2J [T3J [TSJISJ HI |3], involutive 
(n— l)-planes are at the very heart of the theory of characteristics, since the former 
are, in a sense, dual to the latter. Such a relationship between singularity equations 
and characteristics, firstly observed by Krishchenko [7J [5] and Hermann, admits 
a transparent and simple formulation (Corollary |21[ ) in terms of involutive flags. 
These are an obvious generalization of the equation of involutive planes (Definition 



14 1 , and a special cases of integral flags, already appearing in the theory of exterior 
differential systems (e.g., in the proof of the Cartan-Kahler theorem |2J). 

Grassmann and flag bundles (Section [TJ have been put at the basis of our analysis 
not for merely aesthetic reasons. Indeed, this allows to construct a polar distribution 
on Xfc(C) by standard manipulations of multi-linear algebra (Section |2j. We called 
it polar since, quoting Bryant&Grifhths |2j, is related to the linearization process 
in multi-linear algebra known as "polarization", the most common example being 
the polarization of a quadratic form on a vector space to produce a bilinear form. 
Examples of polar distributions (D, in this note) can be found in the context of 
Mongc- Ampere equations (to the authors knowledge, firstly by Vinogradov [41 ITU 
HIJ) and the idea behind it already appears in Bryant&Griffiths' classical book (2 j 
(in the guises of polar space or space of integral enlargements), but it has never 
been studied in full generality. This task, which in fact poses no serious threat, is 
carried out here, in the framework of Grassmann bundles. The subtle point is the 
natural character of the construction, which, though never stressed throughout the 
paper, is its most important feature: it implies, in particular, that D is a contact 
invariant of the nonlinear PDE it is associated with. 

Besides Corollary [2TJ the main results of this note are contained in the last 
Section [5j where it is proved, in some cases, that D is isomorphic to the Cartan 
distribution on a jet space with less independent variables. This fact allow, for 
instance, to associate with a fc th order PDE in 2 independent variables, a family of 
ODEs which is, by construction, a contact invariant of the given PDE. 



Notations and conventions. All vector spaces are over R. If P is a C 00 (M)-modulc, 
we replace C°°(M) by M in all statements involving P (e.g., duals, tensor products, 
homomorphisms) . Also P x denotes the fiber of the bundle corresponding to P, at 
the point x £ M. The same notation is used for non-linear bundles: e.g., E x is the 
fiber of 7T : E ->• M at x £ M . If F : N ->• M is a smooth map, F°(tt) denotes the 
induced bundle on M. The dual of vector spaces (or module) P is denoted by P v . 
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All manifolds and maps are assumed to be smooth and all PDEs are supposed 
to be formally integrable ones. Symbol FT denote the tangential map associated 
with F. We adopt Einstein notation for repeated indexes. 

Concerning definitions and main structures of jet spaces and nonlinear PDEs, 
we adhere to the notation of |T|. 

1. Preliminary results about Grassmannian bundles 

If V is an (n+m)-dimensional vector space, Gr ( V, n) denotes the nra-dimensional 
smooth manifold made of n-dimensional linear subspaces of V. The short exact 
sequence of vector bundles 

(1.1) R{V,nf >Gr(V,n) X V N(V,n) 





Gr(V,n), 

where ttq^ is the trivial bundle, R is the (n-dimensional) tautological (or canonical) 
bundle, and ./V is the (m-dimensional) normal bundle, is known as the universal 
sequence over Gr (V,n). By definition, Rl = L (hence the name "tautological") and 
Nl = 77, for all L e Gr(V,n). Together, the tautological and the normal bundle 
allow to describe the tangent bundle t Q r of Gr (V, n), since 

T Gr = Horn (R, N) = i? v ® JV 

canonically. In order to introduce the characteristics of a distribution on Gr (V, n), 
it is convenient to prove the next Corollary [T] 

Corollary 1. A projective sub-bundle T Gr exists, such that P(t Gt ) is the image of 
the Veronese embedding 

HR y ) XGr(v,n) P(N) CP(r Gr ). 

Proof. Just observe that, thanks to ([TJ, one can speak of the rank of a tangent 
vector to Gr (V, n), which is the rank of the corresponding homomorphism from R 
to N. Then 

71 Gr (V, n) = f {6 E T L Gr {V, n) \ rank (0) = 1} 

is a projective variety in T^Gr (V, n) for all Lg Gr (V, n). Define Tq t as the bundle 
whose fiber at L is (JT|). □ 

Let now A be a distribution on Gr (V, n) , regarded as a vector sub-bundle A C 
TQ r • Hence, P(A) is a smooth sub-bundle of P(rQ r ), and it makes sense to consider 
the sub-bundle 

(1.2) P(4)nP(A) 

of P(A) , made of rank-one directions belonging to A. Denote by pi : P(i? v ) x Gr i Vn \ 



¥(N) — > P(i? v ) the canonical projection, and take the image of (1.2) via p\. The 
so-obtained projective sub-bundle 

char R (A) " Pl (P(4)nP(A)) 

of P(i? v ) is known as the bundle of characteristics of A. 

In order to pass to the global point of view, one just needs to "replace" vector 
spaces by vector bundles. To make this rigorous, it is sufficient to observe that the 
universal sequence ( |1.1[ ) is well-behaved w.r.t. linear transformations of V, i.e., each 
4> € GL (V) induces a diffeomorphism <fi of Gr (V, n), and bundle automorphisms of 
R(V,n), Gr (V, n) x V, and N(V, n), which cover <j>. 
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So, to any vector bundle £ : E — ► M one can associate a smooth bundle Gr (E, n) 
over M , and a short exact sequence of vector bundles over Gr (E, n), 



(1.3) 



R(E, n)<- 



Gr (£, n) X M E ■ 



N(E,n) 



Gr (E, n) 



such that Gr (E,n) x = Gr (E x ,n) and the restriction of (1.3) to a point x G M 
equals ( |1.1| ), with V = E x . Furthermore, if VGv (E, n) denotes the vertical tangent 
bundle of Gr (E, n), then 

VGr (E,n)= Horn Gr {E , n) (R(E, n), N(E, n)) = i? v (E, n) ® Gr (B , n) N(E, n) . 

Finally, for any vertical distribution A on Gr (E, n), i.e., a vector sub-bundle A < 
VGv (E, n), a projective sub-bundle 



(1.4) 



char R (A) < F(R(E,n) v ) 



exists, such that char R (A)| Gr ( B 



char R (A| Gr(£ 



It is natural to call Gr (E, n) the Grassmann bundle (of £), and ( 1.3 1 its univer- 
sal sequence. Accordingly, char A is the bundle of characteristics of the vertical 
distribution A. 

Example 2 (Jet bundles). Let £ = tm be the t ange nt bundle to M, J 1 (M, n) :— 
Gr (TM, n) and TM := tt°_ (TM), and rewrite (1.3| by decorating R and N with 
an index "1". The so-obtained diagram 



R 1 (M, n 




TM 



■N\M, n) 




J l {M, n) 



M. 

contains, besides the definition of the 1 st order jet bundle of n-dimensional 
submanifolds of M, also the main structures necessary to deal with 1 st order 
nonlinear PDEs and to define higher-order jet spaces as well. In the theory of 
PDEs, the canonical and the normal bundles R 1 (M,n) and N 1 (M 1 n) are some- 
times called "distributions" (like, e.g., in pQ), and we embrace this terminology. 
Indeed, R 1 (M, n) can be interpreted as a relative distribution, in the sense that it 
is generated by relative vector fields (see [T5]) w.r.t. 7r 10 (hence, U R" stands for "rel- 
ative"). Also A rl (M, n) can be seen as a distribution, since, for any n-dimensional 
submanifold L C M, the restricted distribution A^ 1 (M, is isomorphic to the 

normal (in the usual sense) distribution of L. 

In this perspective, the Cartan distribution C 1 plays the role of the "absolute" 
counterpart of the relative distribution R (M, n), i.e., 



(1.5) Cl = (d8ir lfi )-\R 1 (M,n) e ), 9 G j\M,n). 

Example 3 (Characteristics of a PDE). When £ is a l st order PDE, i.e., a closed 
sub-bundle £ C J 1 (M, n), then the vertical distribution V£ over £ can be seen as 



the restriction to £ of a vertical distribution A on J X (M, n). Then, thanks to (1.4) 
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above one can speak of the bundle of characteristics of E Q unequivocally defined 
as char R (£ ) := char R (A)|£. 

2. The polar distribution on isotropic planes 

We collect here some necessary algebraic facts about bilinear forms. V and 
W are finite-dimensional vector spaces and, given a subspace L C V, U is the 
annihilator of L in Horn (V, W), i.e., the linear subspace of Horn (V, W) given by 
those homomorphisms which vanish when restricted to L. 

Let n g V v A V v <g) W be a VF-valued skew-symmetric bilinear form on V, take 
an element L G Gr(V,n), and fix an its complement L c . Then n determines a 
"restriction to the graph of h" map 

Horn (L, L c ) L v AL V ®W, 

h — > no(id L + /i), 

which is a polynomial of degree 2 on Horn (L, L c ) with values in L v A L v <E)W since, 
Vx,y G L, 

(e u (h))(x, y) = U(x, y) + Tl(x, h(y)) + U(h(x),y) + U(h(x), h(y)). 
Hence, the differential <ioen coincides with the linear part of |2]), viz. 
(2.1) d e n : h E Horn (L, L c ) i — >Il(-,h(-))+ U(h( ■ ), • ) € i v A L v <g) W. 
Now observe that the subset 

z„(n) d = {Le&^n) I n| L = 0}, 

made of (II-) isotropic n-planes in V, is a quadric variety, since it coincides with 
e r[ 1 ({0}) in the affine neighborhood Horn (L, L c ). In other words, en can be seen as 
a smooth map on Gr (V, n), which is defined only on Horn (L, L c ), where it becomes 
polynomial of degree 2, and its differential at L is given by ( |2.1| . Interestingly, if L 



is isotropic , an analogous formula as (2.1) for d^eri can be given, which does not 
require choosing a complement of L. To this end, for any L G Gr denote by 
the homomorphism 

L — -> V v W, 
I U*(l):=R(l,-), 

and notice that L G I ra (n) if and only if imll* C = (x") V ® W- Hence, if L is 
isotropic, then d^eri reads 

T L Gr {V, n) = Horn (l, j\ I v A L v ® W, 

(2.2) ft ((i >y )—^n#(a:)(M»))-I^(w)(ft(a:))). 

For any L £ I n (IT) put 

T L d = kerd ien = {fee Horn | nf (x)(ft(y)) - nf (y)(ft(x)) = Vx,yG 

D L d ^ {h£Kom R {L,V-)\n#(.)(h(-))=0}. 

Clearly, Dj, is a linear subspace of T^, which in its turn is a linear subspace of 
Ti,Gr (V, n). Hence, T : L i — > Ti, and D : L i — ► £>£ can be seen as distributions 
on Gr (V, n) along the subset I n (Ii). Let I„(n) be the set of smooth points of l n (H), 
i.e., those where (2.2 1 has maximal rank. Then T\f ^ is the tangent bundle of 



I n (n), and D n = D\% , n \ is a distribution on X„(n). 



IX„(H) 



^This is just the "characteristic variety" d la Bryant&Griffiths [2], recast in a bundle— theoretic 
context. 
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Definition 4. D n is the polar distribution^ on l n (H). 

Remark 5 (Geometrical interpretation). Let L be isotropic and h an its deformation, 
i.e, h £ Hom(L, ¥■). This means that h determines a 1-parameter family {L e } of 
n-planes, where 

L e :— graph(e/i) = {x + e/i(x) | a; £ L} 
and Lo = L. Now, condition for L e to be isotropic is 

n| L , = «■ Ol(x, %)) + n(/i(x), y)))e + n(/j(x), %))e 2 = Vz, y G L. 

Hence, h £ Tl if and only if {L e } is to 1 st order made of involutive planes. On 
the other hand, the family {L e } determines an enveloping manifold^] whose tangent 
space at h (seen as a subspace of V) is known as the osculator osc (h) of h, but 
condition h £ is not enough to guarantee that osc (h) is itself isotropic. 

However, h £ Dl if and only if osc (h) has this nice property: for each isotropic 
subspace W C osc (ft), L + W is isotropic as well. This clarifies the geometrical 
content of Definition [4] Dl is made by those deformations of L which are to 1 st 
order orthogonal to L, with respect to n (i.e., it is the "polar space" in the sense of 
0). 

Remark 6 (Duality). Let diml^ = n. In view of the isomorphism Gr(V,n — 1) = 
P(V" V ), the subset I„_i(n) C Gr (V, n — 1) of isotropic hyperplanes in V is in fact 
a projective subspace of P(I^ V ), its equation being a AH — 0, where a £ V v . 

In order to link the purely algebraic considerations of this section with the equa- 
tion of involutive planes which will be discussed below, it is enough to "globalize" 
the above constructions, by starting with a tensor n € T(£ v A £ v ® rf), where 
£ : Q — > M and r\ : Q' — > M are vector bundles. Then X n {U) is a sub-bundle 
of Gr(<5, n) — > Q and D n is a vertical distribution on I n (Ii) In our subsequent 
analysis, Q = C k will be the bundle (over J k (E, n)) of Cartan vectors and n = n fe , 
the mctasymplcctic form on C k . 

3. The equation of involutive planes 

Let A be a distribution on E, and f2 G A 1 (E)® E P, with P a C°°(£')-module, an 
its dual P-valued form, i.e., such that A = kerfi. Denote by n G A 2 (A V ) ®e 
the curvature form of A, where A^ is the 1 st derived distribution. 

Definition 7. A tangent plane R to E is called involutive if 

n\ R = o, n| fi = o. 

The totality of n-dimensional involutive planes of E is the equation of involutive 

n- dimensional planes of E, and denoted by Z„(A). If E is fibered, then I„(A) = 
{R £ I n (A) | R is horizontal} is an open and dense subset of X„(A) called the 
equation of horizontal involutive n-planes. 

A necessary step to reveal that the equation of involutive planes is equipped 
with a canonical distribution, is to show its equivalence with the Grassmannian of 
isotropic planes. 

Lemma 8. I„(A) = I n (Jl) and, in the fibered case, I n (A) C i n (JT). 

Proof. Regard A as a bundle over E, and n as a skew-symmetric bilinear form on 
it, with values in the bundle Q' whose module of sections is (see the end of 
previous Section |2j. Then l n (H) is the sub-bundle of Gr (A, n) made of n-planes 



The name is mutuated by the standard terminology of the theory of quadric surfaces. 
'Called block by some authors. See, e.g., [5]- 
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on which il vanishes or, which is the same, the sub-bundle of Gr (TE, n) = J 1 (E, n) 
made of n-planes on with both £1 and II vanishes, i.e., precisely the equation I n (A). 

It remains to be observed that an horizontal involutive n-plane is always sur- 
rounded by a neighborhood of such planes. □ 

Example 9 (Coordinates). Let E — {(x 1 ,^)}, and A given by means of 1-forms 

A = r~UeAkeru; A , 

with oj a = ujfdx 1 + uj A du j . Then 

(3.1) X n (A)={0G j\E,n) \u; A \ Re =0, du A \ Rlj = 0, \/A e A} 

is locally given by the system of equations {f A = 0, / t 4 = | i, j = 1, 2, . . . , n, Ae 
A}, whereQ 

(3-2) ft = w^+wX 

(3.3) = u^+u^u^+uf^tful 

To prove (3.2) and (3.3) it is enough to recall that a point 9 = (x l ,u a ,uf) corre- 
spond to the plane Rg generated by the n vectors d x i \g + ufd u a \ g, i — 1,2,.. . , n, 
and write down conditions (3.1 ) in terms of the coordinates of 9. 



Remark 10. L et £ C J 1 (-E, n) be given just by (3.2 1 alone. Then 2"„(A) = 
TT2 t i(£ ^'), i.e., (3.3) are differential consequences of (3.2 . 

Example 11. If C k is the Cartan distribution on J k (E, n), then X n (C k ) is the closure 
of J k+1 (E,n) in J 1 (J fe (i?, n), n). Adherent points corresponds to the so-called 
singular i?-planes{^] 



Example 12. Let C k be as in Example[ll| Then 2 n (C k ) = J k+1 (E,n). 

Remark 13. Leaves of 2^°°^ (A) are in one-to -on e correspondence with n-dimensional 
involutive submanifolds of A (see Remark 10). If A is integrable, then Zf° (A) 
can be seen as the orbit space of the leaves of A. 

Lemma[H]above suggests an obvious flag-theoretic generalization of the equation 
of involutive planes, namely 

l n ,„>(A) d = (In (A) x E l n ,(A)) n Gr(A,n,n / ), 
where Gr (A, n, n') is the flag bundle associated with the vector bundle A over E. 

Definition 14. I„ jn /(A) is the bundle of flags of involutive planes of the distribu- 
tion A. 

4. Application to jet bundles 

Definitions and results of Section [3] have interesting applications in the context 
of the geometrical theory of PDEs. Here involutive planes of dimension n — s are 
called of type s. More precisely, the 1 st order PDE 

(4.1) 2„_ s (C fc ) = {0e J 1 (J k (E,n),n-s) | n k \ R(i = 0, U k \ Re = 0} 

(resp., its open and dense subset I„_ s (C fe ) is called the equation of s -type involutive 
subspaces (resp., s-type /lorizoniaZinvolutive subspaces) of C k , where fl k is the TV- 
valued Cartan form and H k its curvature form. 

Corollary 15. Z„_ s (C fe ) is equipped with the natural, i.e., contact-invariant, polar 
distribution D n . 



^Physicists' notations have been used here: "[i, j]" means sew-symmetrization of indexes, and 
; i" means derivative with respect to the coordinate labeled by i. 

Firstly studied by Vinogradov in the context of singular and multivalued solutions 1181 1171 . 
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Example 16. I n (C k ) is the k + 1 st order jet space J k+1 (E,n), and the restricted 
projection 7Ti,olx„(c fe ) ^ s ""fc+i,*- i- e -i i e ^ s of order k + 1 are nothing but involutive 
subspaces of C k of type 0. Notice that, in this case, the polar distribution is trivial. 

Pass now to the next case, i.e., of type 1 subspaces, and introduce the following 
commutative diagram 



(4.2) 



2n,n-l(II *) 



J k (E,n) 



Zn-lOl*- 1 ) 



(E,n) 



Observe that for any 9 e J k (E,n), c e = Gt (R e ,n - 1) = Pi$, i.e, Z„ : „_i(n fe - 1 ) 
is the total space of the bundle Pi? v over J k (E,n). This provides a transparent 
geometrical interpretation of flags of involutive subspaces, and shows that Z„ in _i 



is the natural environment for characteristics (see Corollary 17 below) 



Corollary 17. Let £ be a k order PDE. Then char {£) is a sub-bundle of c\g. 



Diagram (4.2 1 shows the relationship between characteristics and so-called sin- 



gular solutions, since involutive subspaces of type s can be seen as the "shadows" 
of singular i?-planes of type s. Details are as follows. 

Definition 18. Elements of X n {C k ) \I n (C fe ) are called singular R-planes^ An 
n-dimensional, C fe -integral submanifold W of £, is called a singular solution (at 
9 € £)■, iff T S W is a singular i?-plane. The type of a singular i?-plane R < T e J k is 
the codimension in Rg of dgTTk,k-iR- The type of the singularity of W (at 9) is the 
type of T B W. 

Remark 19. The projection Ttk+i,k can be regarded as the tangent map Ttt/c^-i, 
sending horizontal i?-planes on J k into horizontal i?-planes on J k ~ 1 , i.e., acting 
between I n {C k ) and I n (C k ~ 1 ). On the other hand, Tn^.k-i cannot act between 
T n (C k ) and I n (C k ~ 1 ), since the projection of a singular i?-plane is still involutive, 
but fails to be n-dimensional. 

However, if I n (C k ) is seen as a stratified manifold, where the s th stratum is 
given by singular i?-planes of type s, with s = 0, 1, . . . , n, then Tn^.k-i defines a 
projection of the s-stratum of T n (C k ) onto I n -s(C k ~ 1 ) (indeed, the projection of 
an involutive space is always horizontal). 

In other words, 2"„_ s (C fc_1 ), plays the role of "empty equation", i.e., the "recep- 
tacle" for all projections of all singular solutions of type s of all k th order PDEs. 
The next logical step is to confine oneself to a concrete fc th order PDE £, and con- 
sider the elements of I n _ s (C k ~ 1 ) which are the projection of a singular solution of 
£ of type s. More precisely, denoting by £k-i the projection of £ onto J k ~ 1 , the 
projection of the s-stratum of I n (Cg) onto I n ^ s {C^~ 

def 



E W £ = {T7r fc)fc _i(W0 | W € In(Cl) is of type s} C ^(C*^) 



is a 1 st order PDE on £k-i called the s-type singularity equation of £. 

Remark 20. The projection of a regular solution of £ on £k-i, is a solution of 
£k-i, while the projection of a singular solution of £ (of type s at 9) is a singular 
submanifold of £k-i, and the tangent space to its singularity locus at Ttk is 



^This definition is due to Vinogradov 1181 . 
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an (n — s)-dimensional involutive subspace of C £ ~ ± . The totality of such tangent 
spaces is £[ s i£. 

The link between singularity equation and characteristics easily follows from 



Diagram (4.2 1 above 



Corollary 21. is the projection of char £ via the map p of Diagram (4.2 1. 

The next result is a source of contact invariant for nonlinear PDEs. 

Corollary 22. Singularity equation E[ s ]£ is equipped with a natural, i.e., contact 
invariant, polar distribution, obtained by restriction of the polar distribution D n 

5. Computation of D n on some equations of involutive lines 

The equation of involutive lines is particularly simple, since all lines are invo- 
lutive, i.e., Zi(A) = P(A) for any distribution A. Hence, in this case, D n is a 
distribution on the projectivized bundle P(A), where II is the curvature of A. 

While working on the present subject, the authors developed the opinion that, 
if II is the metasymplectic structure on a "empty equation", i.e., a jet space, then 
the corresponding polar distribution D n on it should take a particularly simple 
form. In this context, "particularly simple" means isomorphic to the Cartan distri- 
bution on a PDE of the simplest form, i.e., again an empty one, expectedly with a 
smaller number of independent variables. This guess can be formalized as a simple 
conjecture^ which can be proved true in few cases, collected in this last section. 

Conjecture 1. Let C be the Cartan distribution on J k (E,n), and II its meta- 
symplectic structure. Then an N > n exists such that, for any 9 G J k (E,n), 
the projective space P(Ce) is (locally) isomorphic to J k (R N ,n — 1), and the polar 
distribution D n coincides with the Cartan distribution on J k (M. N ,n — 1). 



In Subsection 5.1 we show that the polar distribution on a 1 st order empty PDE 
in n independent variables identifies with the Cartan distribution on an empty 1 st 
order PDE in n — 1 independent variables. 



In Subsection 5.2 we pass to a nonempty 1 st order PDE £, and show how to 
compute the family of polynomials determining the corresponding polar distribution 
on the equation of the involutive planes of £ , from the smooth function(s) cutting 
out £. This is a nice consequence of Conjecture [T] which can be stated in full 
generality. 

Corollary 23. Let £ C J k (E,n) be a PDE, C £ its Cartan distribution, and 
the corresponding metasymplectic structure. Then, for any 9 € £ , the projective 
space P(Cg ), equipped with the polar distribution D U£ , is a PDE m J k {R N ,n-l). 
Moreover, the so-o btained family {¥(Cl)} ee£ of PDEs is a contact invariant of £ , 
in the sense that any contact transformation F : £ — > £ induces an isomorphism 



between the PDEs P(Cf ) and P(Cp^). Subsection 5.3 outlines a strategy to obtain 
a complete proof of Conjecture [l] to regard k + 1 st order jets as 1 st order PDEs and 
then use the result of Subsection 15.11 as the induction basis and that of Subsection 



5.2 as the induction step. Due to the complexity of the expression involved (see 



(5.141 and (5.15)), the authors did not complete the task, but the subsequent 
examples show it should be conceptually simple. Example [25] proves Conjecture [T] 
for a scalar 2 nd order empty PDE in two independent variables, highlighting the 
role of a good induction step for the general proof. The result is in fact valid 



for all scalar equations in two independent variables (Example 27 1 , but its proof is 



^This conjecture was brought to the attention of the first author by his Ph.D. advisor, prof. 
A. M. Vinogradov. 
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technical, and gives no hints for the general case (for this reason, the proof of its key 
Lemma [28] which is rather lengthy and can be found in the second author's Ph.D. 
thesis, is omitted). In the Example 29 (which is an application of Corollary 23 1, 
results of Example [27] are used to describe the polar distribution on the singularity 
equations of a third order nonlinear PDE in two independent variables, showing 
the importance of Corollary [21] in these kinds of computations. 

5.1. Computation of D n for T\ on 1 st order jets. Let £ be a manifold of 
dimension n + m. Denote by C and II the Cartan distribution on J 1 (i?, n) and the 
metasymplectic structure on C, respectively. Choose local coordinates (x l , , u?) , 
with i G {1, . . . , n} and j G {!,..., m}. Then, locally, 



C = Span{A,d u J 



(5.1) 
and 

LT = dx 1 A du\ ® d u j . 

Since all lines lying in C are II-involutive, 7i(II 1 ) = PC. Hence, D n is a vertical 
distribution on the bundle PC — > J l {E,n). Since J 1 (£',n) is homogeneous w.r.t. 
contact transformations, it is enough to describe the restriction of D n to a generic 
fiber PC e , with 6 G ^{E.n). Let L = Span -{7} G PC e , with 

(5.2) l = D 1 \ e + b a D a \g + fid u] 

* 9 

where a E {2, . . . , n}. Then 

(5.3) (b 2 /,"../V /m 

is a coordinate system for the open neighborhood Horn r (L, L c ) of L, in the projec- 
tive space FCg - 



sn(m+l)- 



, where 
L c := 



Span j-D a |e, j 



is a complement of L. Since Homj{(i, L c ) identifies with T^(PCe) by means of the 
correspondence 



[h:l\-> B a D a \, 

the fact that belongs to 
F( . More precisely, 



■b°-\g 



>■: <>f 



5^ reflects on a condition on the coefficients B a and 



ILg(l,h(l)) = 
& U e (D 1 \ e + b a D 



& F(U[D u d u 



+ ./;;/rn (a , 



= 



+ f>d ul \e,B a D 

b a F^n (p a ,d ui 

(Fl+b a Fi-f a B a )d u i\g=Q 

F{=fiB a -b a Fi Vj = l,2,...,m. 

Equations (5.4 1 shows that the elements of are in one-to-one correspondence 
with the (n— l)m-tuples (B 2 , . . . , B n , F%, . . . , F£) . In particular, we obtain a basis 



(5.4) 



6,-- 

and B a 



; <J| for the X^'s (with /3 = 2, 3, . . . , n), and 



of by choosing F^ 
then choosing B a = and = <Sf<^ for the X%'s (with = 2,3, 
fc = 1, 2, . . . , m). Notice that the values at L of the vector fields 

(5.5) X a := d b . 



and 



(5.6) 



Xf 
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P n is spanned by (5.5 1 and (5.6 1. 



are, by construction, £ Q and respectively, i.e. 
Moreover, 



(5.7) 



[X a ,Xf] 



Observe that (5.5 1 resembles a total derivative, (5.6 1 is almost a vertical vector 
field, and (5.7 1 is reminiscent of the standard commutation relation of the standard 
non-holonomic frame of a first order jet space. The idea is to look for a coordinate 
system which "turns right" formulas (5.5 1, (5.6) and (5.7 1. To this end, let 



(5.8) 


(y 2 


...,y n ,v\ 




v m y 2 ,. 


be a new coordinates on 


PC 


, given by 






(5.9) 




b a := 


y a , 




(5.10) 




fj := 

J a 






(5.11) 




fl ■■= 


2v j 


-<y a 


Then it is immediate to 


see 


that 







v a ^vi — X a , 



[x a ,xf]. 



In other words, (5.9 ), (5.10 1 and (5.111 define a local diffeomorphis between PCg and 
ji n _ l) 5 which sends D 11 into the Cartan distribution of J 1 (M n - 1+m , n- 

1). 

5.2. Computation of D n for T\ on 1 st order PDEs. In the same setting of 
Subsection 5.1 the dual basis of (5.1 ) is 

(5.12) 



C v = Span jtfccSdttA 



Let now £ C ^(E.n) be the PDE determined by e C°°(J l {E,n)), and put 
C £ := C\s. Accordingly, H £ will denote the metasymplectic structure on C £ . In 
terms of coordinates (5.121, := de4>\c e € Cg reads 

a <l> = D i (<p)(6)d eX i + ^ 
ou J - 



(0)d e ul 



Hence, a line I £ P(Cg) belongs to I\ {Cf) = P(ker a^) if and only if a^il) = 0. More 
precisely, if I is given by (5.2 1, then 



(5.13) 



a*(f) = DMW) + b a D a {4>){e) + fA{6) 

dui 



and the right-hand side of ( |5.13[ ) is a polynomial in t he coordinates (5.3 1 which 
determines T\ (Cf ) . Passing to the coordinates (5.8), the main result of Sub- 
section 5.1 guarantees that Zi(Cf) is isomorphic to the equation {P* = 0} in 



jl(K»-TTrn jn _ 1) ) w ith 

Pt{y a yy a ) = Di{<t>)\e + D a m e y a + 2 



- y a vi) + 



du 3 a 



In particular, D n is isomorphic to the restriction of the Cartan distribution on 
J 1 (M n - 1+m ,n- 1). 
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5.2 



can 



5.3. Computation of D n for I\ on higher— order jets. k + 1 st oder jet s can be 
treated as differential equations on k th order jets. So, results of Subsection 
be tested in this context. 

Example 24 (k = 1). Relabel coordinates of E' = J 1 (£ , ,n) as follows 



E' = {x',u J } with J G {!,... m} U ({1, . . .n} x {1, . . .m}) , 



where 



Hence (see Example 12 before), J 2 (E,n) 
where 



{F( = 0, GY = 0} c J 1 ^',™), 



/•;•' := u ; •»/ •'. 

g; j : = uJ* , -«{ J . 

Observe that the Cartan distribution C of J 2 (E, n) is the restriction to J 2 (E, n) 
of that on J 1 (i?', n), and denote by n 2 the corresponding nictasymplectic structure. 
Then (see Subsection 5.2 1, the distribution D n on P(Cg) is isomorphic to the Cartan 
distribution on J 1 (R"~ i+m+nm , n — 1), restricted to the system of equations given 



by the m "(" +1 ) polynomials (the (5.14 1 's are symmetric in i, I) 



(5.14) pf := -u l { j (e)-u^(0)y a + SiS j J (2v J -y a v J a )+SrS j J v J a , 

(5.15) pf := 2(4^-^)^-^) + ^^, 
in the mn(n + 1) +n — 1 variables 



y 2 ,...,y n ,...,v J , 



Example 25 {k = 1, n = 2, m = 1). If g} are local coordinates on 

J 1 (M 3 ,2), then the eleven-dimensional jet space J 1 (J 1 (M 3 , 2), 2) has coordinates 

x,y,u,p,q,u x ,u y ,p x ,p y ,q x ,q y , 

and its eight-dimensional Cartan distribution C is spanned by 

DD8 8 8888 

Hence, the eight-dimensional 2 nd order jet space J 2 (]R 3 , 2) identifies with the equa- 
tion 

(5.16) £ = {F 1 = F 2 = G = 0}, Fx := u x - p,F 2 := u y - q,G := p y - q x . 
Let 

^7 /a; 7 /$/ i Jjc i fy j fx i fy 

be the projective coordinates of the line 

L := Span { D g \ g + b D y \ g + f x x d Ux \ e + ft d Uy \ g + f 2 x 8 Px \ g + f 2 8 Py \ g + f x 8 qx \ g + f y 8 qy \ g } 

in the seven-dimensional projective space P(Cg). It is convenient to present the 
four-dimensional distribution D n by means the three 1-forms 



(5.18) 



u> :=df-Edb + bdE. 



Indeed (see (5.5) and (5.6)), D n is spanned by the four vector fields 

<h, • ./',;'''/ • B fy h <>. ■ 
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which are also a basis for the common kernel of the forms (5.181. In the present 

j) read 



setting, polynomials (5.141 and (5.15) on 

(5.19) p e Fl = fl- Px {e)-bpy{0), 

(5.20) = fi-q x (0)-bqy(0), 



(5.21) 



The zero locus of fl5.19|), (|5.20| and (|5.2l| is precisely P(Cf ) ( as in Subsection 



5.2 



), an d (see Subsectio"n"]5.1[) P(Cg) is isomorphic to J^M 4 , 1). Hence, ( |5.19| , ( |5720p 
and (5.21 1 together define the 1 st order equation 



(5.22) 

(5.23) 
(5.24) 



Lq y (9)+tQ(8) + ~ Px (8), 

Q(9)+tq y (6), 
2v 3 - vft, 



where Q := p y = q x . 

We claim that P(Cf ) is isomorphic to J 2 (M. 2 , 1). To this end, observe that con- 
dition (5.221 is algebraic, i.e., it does not involve derivatives, and that (5.231 is 
an its differential consequence. This means that P(Cf) can be understood as a 
codimension-one submanifold sitting in the five-dimensional graph M of the map 



(t,v 2 ,v 3 ,v 2 ,v 3 ) 



t,jq v (9)+tQ(e) 



- Px (0),v 2 ,v 3 ,Q(9)+tq y (e) 1 v 2 ,v 3 



dv J — v\dt the three 1-forms which define the Cartan distribution C 



Let rf> 

on J 1 (M 4 ,1). Then, by construction, ry 1 ^/ 



0, and, in terms of the coordinates 



(t, v , v , v 2 , vf) on M, C'\m is given by the two 1-forms 



(5.25) rf = dv 2 - v 2 dt 7 jf = dv 3 

In other words, after the change of coordinates 



v 3 dt. 



t 
u 

P 
u t 
Pi 



t. 



M becomes isomorphic to the five-dimensional iterated jet spaces J 1 (J 1 (K 2 , 1) , 1) = 
{(t,u,p,u t ,Pt)}, and C'\m identifies with its Cartan distribution. Since J 2 (R 2 ,1) 
is singled out by the equation u t — p. which, in the new coordinates, reads 

(5.26) vf=v 3 , 



it remains to be pr oved that (5.241 and (5.26) gives the same equation. To this end, 
observe that (5.241 allows to express v 3 as a function of (t, v 2 , v 2 , v 3 ) or, equivalcntly, 
p as a function of (t,u,u t ,Pt), which can then be taken as internal coordinates on 
P(Cg ). Accordingly, after restriction to P(Cf), forms (5.251 become 



(5.27) 
(5.28) 



m 



dti 
1 
2 



u t dt, 
d(u t +Ptt) 



Ptdt = \{du t - Pt dt + tdp t 



The fact that the distribution determined by (5.271 and (5.281 is the same as the 
Cartan distribution on J 2 (M 2 , 1) is a consequence of the general Lemma 26 below. 
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Lemma 26. Let A be the distribution on M 4 = {(t, u,p, q)} given by the two forms 

rji = du — pdt, 

772 = dp — qdt — atdq, 

with a G R\{1, i}. Then A is isomorphic to the Cartan distribution on J 2 (R 2 , 1). 

Proof. With the new coordinates 

t = t, 

u = (1 — 2a)u + apt, 
p = atq + (1 - a)p, 

q = q, 



one obtains 



and 



rji = (1 — 2a)du + ad{pt) — {atq + (1 — a)p)dt 
= (1 — 2a)(du — pdt) + at(dp — qdt) 

V2 — d(atq + (1 — a)p) — qdt — atdq 
= (1- a) (dp -qdt). 



a 



Example 27 (k > 1, n = 2, m = 1). Let (x, y, u a ^ \ a + b < k) be lo cal coordinates 
on J fc (M 3 , 2), and C its Cartan distribution. Much as (5.3 1 and (5.17), the projective 
coordinates 

(5.29) 6,/ ,/i,...,/fc 

on the (k + 2)-dimensional projective space P(Cg) identify the line generated by 

k 

l-.= D x \ e + bD y \ e + Y,h d ^ h -i\e- 

1=0 

The metasymplectic form on C reads 

k-l 

n = ^2(dx A dui+i ® d Ul k _ l +dyA du k -i,i+i ® d Uk _ u i) 
2=0 



and (see (5.18 1) the distribution D n is the common kernel of the k one-forms 
(5.30) Ul = df l+1 - fidb + bdf h 1 = 0,..., k-l. 



By Lemma 28 below, after a suitable change of coordinates on V(Cg), the distribu- 
tion D n coincides with the Cartan distribution on J k (M. 2 , 1) = {(y, v, V 1 , . . . , v k )}. 

Lemma 28. Change of coordinates 

b = y 



, „^ (j + i)H/)'-j , , 

h = 11 1^ (T—)l v i l = °i--- k 

j=o v 



makes the system (5.301 equivalent to 

dv-' — v^ +1 dy. 
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Example 29. The equation £ — {(/) — 0}, with 

(f> := Uyyy ~ U xxy + U XXX U X yy, 



is a sub-bundle of 

J 3 := J 3 (R 3 ,2) — >J 2 (R 3 ,2) =: J 2 . 

Indeed, </> does not depend on the coordinates of J 2 and as such, for any 8 € J 2 , it 
can be regarded as a function on the fiber J| which singles out £§, viz. 

(5.31) £ s = {cf> = 0} C J$. 

Notice that £g is a quadric hypersurface in the four-dimensional space J|. In the 
sequel, we denote by 8 the projection of 8 on M, and by 8 6 J 3 any point of Jf. 
Let now i7 := Span{w} g P(i?~) be the pencil determined by a covector g 

By definitional £1 is a characteristic of £ at if and only if 

(5.32) 3e€ Afcs.t. u®e€ 



It is well-known^] that condition (5.32) above is equivalent to the fact that the third 



symmetric power or, which belongs to S (Tg M) = TJjj, actually sits in Tg£g. In 
other words, 

charf (£) = {n e P(i$) \ uo 3 € Tj£ e }. 
Notice that F(R~) is one-dimensional, so that correspondence 

6 

(5.33) a f> Sl„ := Span {— ad-gx + dg-y} 

is enough to define a coordinate system on P(i?~). Accordingly, to each value 
of a it corresponds a vertical tangent vector £ a , namely the one associated with 
(—adgx + dgy) 3 via the identification S 3 (TgM) — TgJg, 

f a := — a 3 d u k + 3a 2 3 u \-x— 3a d u + d u 

Since £ a belongs to Tg£g if and only if £ a (0) = 0, the subset char~(£) C P(i?~) can 
be (locally) identified with the set of the zeros of the polynomial 

Pg(a) := u xy y(8)a 3 - 6u xxy (8)a 2 + 3u xxx (8)a - 1. 



21 



Passing to the the fold-type singularity equation £m£ C J 1 ( J 2 , 1), Corollary 
shows that 

(5.34) (E[i]£)e = {kerfJ a | pg-(a) = for some 8 e £ e } C P(C e ), V6» € J 2 , 

where by the kernel of a pencil of covectors we mean that of any its representative. 
Observe that ker tt a is an element of the projective line P(ifo), canonically identified 
with a subset of the four-dimensional projective space V(Cg). In other words, 
(£m£ )e is the union of all the "dual lines" to the characteristics of £ at all points 
projecting over 8, i.e., ( |5.34[ ) can be rewritten as 

(5.35) U U kerr> a C |J P(i? ? ) CP(C e ), G J 2 , 

which seems to suggest that, in order to compute the singularity equation at 8 one 
has to find all the roots of the polynomials pg, with 8 G £g. In fact, in the present 
example, this step can be circumvented. 



3 See |IJ and Example [i] 
'See, for instance, [2]. 
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In view of the definition (5.331 of Q a , its "dual line" is the element of P(Rg) 
represented by r\ a :— D x \g + aD y \g and ( |5.35 ) reads 

(? [1] S)o= IJ U S P an {^}, We J 2 . 



(5.36) 



Equipping P(Cg) with the same projective coordinates (b, fo, /i, fz) as Example 27 
above, Span {rj a } is given by 



b 
fo 
fx 
h 

Observe that A := u xxx {0), [i 



a, 

{6) + bu yyy (9), 
(9) + bu xyv (9), 
{0) + bu xxy (9). 

1 xxy(&) and v '.— u xyy (9) can be taken, thanks 



to (5.31), as coordinates on Eg, and as such their values can be arbitrary. Hence, 
(Em£)e is the three-dimensional projective subvariety of P(Cg) given, in the do- 
main of definition of the coordinates (b, fo, f\, f 2 ), by the mixed implicit/parametric 
equations 

z/& 3 - 6m& 2 + 3A6 - 1, 
v + &(/i 2 — 
jtx + bv, 
A + 



(5.37) 
(5.38) 
(5.39) 
(5.40) 





fo 

h 
h 



where A,/i, v € K. Observe that the first equation corresponds to the innermost 
union appearing in (5.361, while the last three equations correspond to the outer- 
most one. As announced, system above can be brought in a purely implicit form. 
Indeed, (5.381, (5.391 and ( 5.40[ ) allow to express A,/i, v in terms of fo,fi,f2, 

(h + hf2)b-h 
M fib 2 + f 2 b-l ' 

f?b-f 



A = h- 



fib 2 + f 2 b-V 

((/o + /i/ 2 )b-/i)6 
/i& 2 + hb - 1 ■ 



These expressions, substituted in (5.371, give the polynomial equation 
1 + 8/16 2 + flb A + (3/ 2 6 - Qhb 2 - 4)/ 2 6 - 10/ & 3 = 



of the three-dimensional projective variety (EmfQg C ¥(Cg). Finally, performing 
the change of coordinates of Lemma 28 above, we get the 2 nd order ODE 

l+10/u 2 +24y 4 i}+108y 8 ii 2 +4(8-16y 3 w+25y 2 «+54y 4 i;)y 2 u-2/o(ll+36y 4 w)y 3 



0, 



which is by construction contact invariant of £ . 



Acknowledgements. It is a pleasure for the second author to thank to the the Grant 
Agency of the Czech Republic (GA OR) for financial support under the project 
P201/12/G028. 



References 

[1] A. V. Bocharov. . . [et al.], I. S. Krasil'shchik (editor), and A. M. Vinogradov (editor). Sym- 
metries and conservation laws for differential equations of mathematical physics, volume 182 
of Translations of Mathematical Monographs. American Mathematical Society, Providence, 
RI, 1999. 



ON THE POLAR DISTRIBUTION FOR SINGULARITIES EQUATIONS OF NONLINEAR PDFJS 7 



[2] R. L. Bryant, S. S. Chern, R. B. Gardner, H. L. Goldschmidt, and P. A. Griffiths. Exterior 
differential systems, volume 18 of Mathematical Sciences Research Institute Publications. 
Springer- Verlag, New York, 1991. 

[3] D. Catalano Ferraioli and A. M. Vinogradov. Differential invariants of generic parabolic 
monge-ampere equations. Diffiety Preprint Series, DIPS-1/2008, 2008. 

[4] Alessandro De Paris and Alexandre M. Vinogradov. Scalar differential invariants of symplectic 
Monge-Ampere equations. Cent. Eur. J. Math., 9(4):731-751, 2011. 

[5] I. S. KrasiPshchik, V. V. Lychagin, and A. M. Vinogradov. Geometry of jet spaces and non- 
linear partial differential equations, volume 1 of Advanced Studies in Contemporary Mathe- 
matics. Gordon and Breach Science Publishers, New York, 1986. 

[6] A. P. Krishchenko. The structure of the singularities of the solutions of quasilinear equations. 
Uspehi Mat. Nauk, 31(3 (189)):219-220, 1976. 

[7] A. P. Krishchenko. Foldings of JJ-manifolds. Vestnik Moskov. Univ. Ser. I Mat. Meh., (1):17- 
20, 1977. 

[8] F. Lizzi, G. Marmo, G. Sparano, and A. M. Vinogradov. Eikonal type equations for geomet- 
rical singularities of solutions in field theory. J. Geom. Phys., 14(3):211-235, 1994. 
[9] V. V. Lychagin. On singularities of the solutions of differential equations. Dokl. Akad. Nauk 
SSSR, 251(4):794-799, 1980. 

[10] V. V. Lychagin. Geometric singularities of the solutions of nonlinear differential equations. 
Dokl. Akad. Nauk SSSR, 261(6):1299-1303, 1981. 

[11] V. V. Lychagin. Characteristic classes of solutions of differential equations. Dokl. Akad. Nauk 
SSSR, 271(6):1320-1324, 1983. 

[12] V. V. Lychagin. Singularities of multivalued solutions of nonlinear differential equations, and 
nonlinear phenomena. Acta Appl. Math., 3(2):135-173, 1985. 

[13] V. V. Lychagin. Geometric theory of singularities of solutions of nonlinear differential equa- 
tions. In Problems in geometry, Vol. 20 (Russian), Itogi Nauki i Tekhniki, pages 207-247. 
Akad. Nauk SSSR Vsesoyuz. Inst. Nauchn. i Tekhn. Inform., Moscow, 1988. Translated in J. 
Soviet Math. 51 (1990), no. 6, 2735-2757. 

[14] Michal Marvan, Alexandre M. Vinogradov, and Valery A. Yumaguzhin. Differential invari- 
ants of generic hyperbolic monge-ampere equations. Cent. Eur. J. Math., 5(1):105-133 (elec- 
tronic), 2007. 

[15] Jet Nestruev. Smooth manifolds and observables, volume 220 of Graduate Texts in Mathe- 
matics. Springer- Verlag, New York, 2003. 

[16] D. J. Saunders. The geometry of jet bundles, volume 142 of London Mathematical Society 
Lecture Note Series. Cambridge University Press, Cambridge, 1989. 

[17] A. M. Vinogradov. Many-valued solutions, and a principle for the classification of nonlinear 
differential equations. Dokl. Akad. Nauk SSSR, 210:11-14, 1973. 

[18] A. M. Vinogradov. Geometric singularities of solutions of nonlinear partial differential equa- 
tions. In Differential geometry and its applications (Brno, 1986), volume 27 of Math. Appl. 
(East European Ser.), pages 359-379. Reidel, Dordrecht, 1987. 

[19] A. M. Vinogradov. Symmetries and conservation laws of partial differential equations: basic 
notions and results. Acta Appl. Math., 15(l-2):3-21, 1989. Symmetries of partial differential 
equations, Part I. 

[20] A. M. Vinogradov. Cohomological analysis of partial differential equations and secondary 
calculus, volume 204 of Translations of Mathematical Monographs. American Mathematical 
Society, Providence, RI, 2001. 

[21] A. M. Vinogradov, M. Marvan, and V. A. Yumaguzhin. Differential invariants of hyperbolic 
Monge-Ampere equations of general form. Dokl. Akad. Nauk, 405(3) :299-301, 2005. 

Mathematical Institute in Opava, Silesian University in Opava, Na Rybnicku 626/1, 
746 01 Opava, Czech Republic. 

E-mail address, G. Moreno: Giovanni.MorenoOmath.slu.cz 

Lucerne University of Applied Sciences and Arts,, Werftestrasse 4,, CH-6002 
Luzern, Switzerland. 

E-mail address: michael .baechtoldShslu. ch 



